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Ising spin glasses with bimodal and Gaussian near-neighbor interaction distributions are studied 
through numerical simulations. The non-self-averaging (normalized inter-sample variance) parame¬ 
ter U22{T,L) for the spin glass susceptibility (and for higher moments Unn{T, L)) is reported for di¬ 
mensions 2, 3, 4, 5 and 7. In each dimension d the non-self-averaging parameters in the paramagnetic 
regime vary with the sample size L and the correlation length ^(T, L) as UnniP, L) = [Kd£,{T, L)/L]‘^, 
and so follow a renormalization group law due to Aharony and Harris [H]. Empirically, it is found 
that the Kd values are independent of d to within the statistics. The maximum values [Unn{T, Z/)]max 
are almost independent of L in each dimension, and remarkably the estimated thermodynamic limit 
critical [[/„„(T, L)]niax peak values are also dimension-independent to within the statistics and so are 
’’hyperuniversal”. These results show that the form of the spin-spin correlation function distribution 
at criticality in the large L limit is independent of dimension within the ISG family. Inspection of 
published non-self-averaging data for 3D Heisenberg and XY spin glasses the light of the Ising spin 
glass non-self-averaging results show behavior incompatible with a spin-driven ordering scenario, 
but compatible with that expected on a chiral-driven ordering interpretation. 


I. INTRODUCTION 

The non-self-averaging parameter, usually noted A 
or [/22, represents the normalized inter-sample variabil¬ 
ity for systems such as diluted ferromagnets or spin 
glasses where the microscopic structures of the interac¬ 
tions within individual samples are not identical. The 
parameter is defined for ferromagnets as the inter-sample 
variance of the susceptibility normalized by the mean sus¬ 
ceptibility squared [l|, 


C^22(/3, L) 


( L) y ^ var((g")) 

\xiT,L)J [{q 2 )f 


( 1 ) 


where a{T,L) is the standard deviation of the equilib¬ 
rium sample-by-sample distribution of the susceptibility. 
We denote by (• • •) the thermal mean for a single sample 
and by [• • •] the sample mean. In ISGs the spin glass sus¬ 
ceptibility replaces y. The non-self-averaging definition 
can be widened to other observables [I| ; we will also dis¬ 
cuss the behavior of non-self-averging of higher moments 
and (g^) of of the spin-spin correlation q. 

Aharony and Harris [I[ gave a fundamental renormal¬ 
ization group discussion of non-self-averaging in diluted 
ferromagnets, which can be applied also to spin glass 
models. First, they showed that in the paramagnetic 
regime, at temperatures above the critical temperature, 
U 22 (which they referred to as R^) behaves as 


U22{T,L)r^{P{T,L)/Lf (2) 


where d is the dimension of the system. This rule can 
be understood on a simple physical picture : the inter¬ 
sample variability depends on the ratio of the sample 
volume to the correlated volume. Roughly, each sample 


is contained in a ’’box” of volume L'^. When this box vol¬ 
ume is much larger than the correlated volume ^(T)^, all 
samples will have essentially identical properties; when 
the inverse is true, each sample has its own individual 
properties. 

Then at the critical point where £,{T) diverges in 
the thermodynamic limit ThL, 1/22(Tc,L) becomes in¬ 
dependent of L even when L tends to infinity [l[. In 
this strongly non-self-averaging regime the observables 
for each individual sample have different properties. The 
passage as a function of temperature in the thermody¬ 
namic limit from ’’all samples identical” (randomness 
irrelevant) to ’’all samples different” (randomness rele¬ 
vant) is a fundamental signature of the physical mean¬ 
ing of ordering in systems with disorder or in spin-glass¬ 
like systems. Aharony and Harris show that the value of 
U 22 {Tc, L) in the limit of large L should be universal, for 
ferromagnets with different forms of disorder in a given 
dimension. We find empirically that within the ISG fam¬ 
ily this critical parameter is dimension-independent, i.e. 
hyperuniversal. 

We report non-self-averaging measurements in near 
neighbor interaction Ising spin glasses (ISGs) having di¬ 
mensions 2, 3,4, 5 and 7, with bimodal or Gaussian near 
neighbor interaction distributions. There is considerable 
regularity in behavior throughout all this range of d, 
which includes the special cases d = 2 where = 0, and 
d = 7 which is above the upper critical dimension d = 6. 
In the paramagnetic regime U 22 {T,L) = [Kd£,{T,L)/LY 
with Kd ~ 2.5 for all d studied, to within the statistical 
accuracy. Secondly, the peak in 1/22 (T,L) as a function 
of T for fixed L has a value 1/22 (max) for each L which, 
after weak small size effects, is independent of L and also 
of d to within the statistics, ?722(max) « 0.205. The lo¬ 
cation of the peak T(1722(max)) approaches Tc from the 
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paramagnetic regime (higher T) for d < 4 and from the 
ordered regime (lower T) for d > 4. The same rules are 
followed for the higher moments of the spin-spin correla¬ 
tions. 


II. SIMULATIONS 

The standard ISG Hamiltonian is 

JijSiSj (3) 

ij 

with the near neighbor symmetric distributions normal¬ 
ized to {Jij) = 1. The normalized inverse temperature is 
/3 = The Ising spins live on simple hyper- 

cubic lattices with periodic boundary conditions. The 
spin overlap parameter is defined as usual by 

i 

where the sum is taken over all spins and A, B indicate 
two copies of the same system. The spin glass suscepti¬ 
bility is then defined as usual x(/3, L) = 

The equilibration techniques (which are different in di¬ 
mension 2) are described in Refs. Si- On a technical 
level, it turns out that the values of U 22 and particu¬ 
larly the peak value can fluctuate slightly in an irregular 
manner as at each size they depend sensitively on strict 
equilibration having been achieved. This can be used as 
a convenient empirical test for equilibration. 


III. DIMENSION 2 

It is well established that short range ISGs in dimen¬ 
sion 2 only order at T = 0 [1, Q. The Gaussian ISG 
has a non-degenerate ground state and a continuous en¬ 
ergy level distribution. The bimodal ISG has an effec¬ 
tively continuous energy level regime down to an L de¬ 
pendent cross-over temperature T*{L) below which the 
thermodynamics are dominated by the massively degen¬ 
erate ground state This is a finite size regime; in 
the thermodynamic limit regime the bimodal ISG can be 
considered to have an effectively continuous energy level 
distribution similar to that of the Gaussian ISG. 

Measurements on two bimodal models and the Gaus¬ 
sian model ISG in dimension 2 Q show a clear scaling of 
U22{T, L) as a function of ^(T, L)/L, with all the maxima 
in U 22 {T,L) close to 0.20. We show for the standard bi¬ 
modal ISG in dimension 2, Fig.[Tl the data scaled against 
^(T, L)/L on a log-1^ plot. This brings out the fact (not 
mentioned in Ref. [Tj) that for temperatures above the 
peak location temperature, the Aharony-Harris rule [l| 
U22{T,L) = {K2£,{T,L)/Lf holds, with K2 = 2.5(1). 
Below the peak obvious finite size effects due to the 
crossover to the ground state dominated regime set in. 



FIG. 1. (Color on line) Bimodal 2D ISG. Non-self-averaging 
parameter U 22 {T, L) against the normalized correlation length 
\{T,L)/L. L = 12, 16, 24, 32, 48, 64, 96, 128 from right to 
left. The straight line has slope 2. 


From the same simulation runs, data for the higher 
moments {T,L) and {T,L) were obtained and 
the values of the normalized variances U^-i(T,L) and 
U 44 {T,L) were evaluated. Equivalent plots to Fig. [T] 
are shown for U 33 {T,L) and U 44 ,{T,L) in Figs. [D and 
[3] with U 33 {T,L) = (3.29C(T,L)/L)2 and U44{T,L) = 
(4.36^(r,L)/L)2. 

The same data are presented as 1/22 (T,L) against T 
for fixed L in Fig. 01 the peak location is moving towards 
T = 0 with increasing L, and the maximum value is very 
gradually growing with increasing L. A simple extrapo¬ 
lation of the peak data from L = 4 to L = 128 indicates 
a limiting infinite L peak value close to 0.200. 

The U 33 {T,L) and U 44 {T,L) peak values evolve in a 
very similar way to the 1/22 (T,L) peaks, extrapolating 
to large L limit values U33 = 0.38(1) and U44 = 0.60(1), 
Figs. [5] and El On the low temperature side of the bi¬ 
modal data, a minimum in each of the Unn (F, L) at an 
L dependent temperature followed by a plateau (see 0) 
provides a clear indication of the crossover from the ef¬ 
fectively continuous energy level regime to the degener¬ 
ate ground state dominated regime. For the largest sizes, 
this crossover lies below the lowest temperatures at which 
measurements were carried out. 

Data for 1/22 (F,L) for the Gaussian ISG (not shown) 
are very similar to the bimodal data, except that there 
is of course no crossover effect. 

The zero temperature infinite size limit can be defined 
in two ways. Taking the successive limits L —>■ 00 , T 
0 gives an extrapolated value 1 / 22 ( 0 , 00 ) = 0.200(5) for 
both bimodal and Gaussian models, while the successive 
limits T —>■ 0, L ^ 00 gives a value ss 0 in the Gaussian 
case; on the present data it is hard to estimate in the 
bimodal model. 
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FIG. 2. (Color on line) Bimodal 2D ISG. Non-self-averaging 
parameter U33{T, L) against the normalized correlation length 
i{T,L)/L. L = 24, 32, 48, 64, 96, 128 from right to left. The 
straight line has slope 2. 


FIG. 4. (Color on line) Bimodal 2D ISG. Non-self-averaging 
parameter 1/22 (T, L) against the temperature T. L = 128, 96, 
64, 48, 32, 24, 16, 12, 8, 6, 4 from left to right. The straight 
line extrapolates to criticality at T = 0. 




FIG. 3. (Color on line)Bimodal 2D ISG. Non-self-averaging 
parameter Uii{T, L) against the normalized correlation length 
^(T,L)/L. L — 24, 48, 64, 128 from right to left. The straight 
line has slope 2. 


IV. DIMENSION 3 

The bimodal ISG in dimension 3 has a transition tem¬ 
perature for which the most recent estimate is = 
1.102(3) [i^ , and the Gaussian ISG has a transition 
temperature estimated to be = 0.951(9) Q. The crit¬ 
ical values of the dimensionless correlation length ratio 
[^{T,L)/L\c are estimated to be 0.652(3) and 0.635(10) 
respectively. 

Hasenbusch, Pellisetto and Vicari Q have generously 
posted their raw tabulated simulation data for the bi¬ 
modal ISG in dimension 3 on the EPAPS site correspond¬ 
ing to their publication. In addition to the present mea¬ 
surements on 2^^ samples of sizes L = 4,6,8,10,12 we 
have extracted a selection of values of t/22(/3, L) from the 


FIG. 5. (Color on line) Bimodal 2D ISG. Non-self-averaging 
parameter 1/33 (T, L) against the temperature T. L = 128, 96, 
64, 48, 32, 24, 16, 12, 8 , 6 from left to right. The straight line 
extrapolates to criticality at T = 0. 


tables of @ , choosing the data sets with the largest num¬ 
bers of temperatures, L = 16, 20, 24. In each case the 
data correspond to measurements on about 10^ samples. 

The U22{T,L) bimodal data in 3D have almost L- 
independent peak values t/ 22 (T, L)max = 0.207(3) with 
peak locations tending gradually downwards towards Tc 
as L increases. Fig. |8] (see who observed also a very 
similar peak height for a next-nearest-neighbor model). 
Small fluctuations as a function of L can be put down to 
residual equilibration differences as the statistical errors 
in these data are very small because of the large num¬ 
bers of samples. At the critical temperature the finite 
size scaling limit is U 22 {Tc,L) = 0.147 

When scaled against ^{T,L)/L, in the paramag¬ 
netic range 1 / 22 ( 7 ", A) = [K3^{T, L)/L]^ following the 
Aharony-Harris law, with K 3 = 2.6(1), Fig. [3 The peak 
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FIG. 6. (Color on line) Bimodal 2D ISG. Non-self-averaging 
parameter UiA{T,L) against the temperature T. L = 8, 16, 
32, 64, 128 from right to left. The straight line extrapolates 
to criticality at T = 0. 



FIG. 7. (Color on line) Bimodal 3D ISG. Non-self-averaging 
parameter 1/22 (T, L) against the normalized correlation length 
^{T, L)/L. L = 10, 12, 16, 20, 24 (pink, purple, cyan, green, 
black) right to left. L = 16, 20, 24 from [^. The straight line 
has slope 3. 


locations correspond to ^{T,L)/L Ki 0.35. 

In the large L limit, the C /33 and C /44 peak locations 
are tending to Tc, and the peak values extrapolate to 
C /33 ~ 0.39 and C /44 0.61, Fig. 1^ and [TUI 


V. DIMENSION 4 

C/ 22 (F, L), C/ 33 (r, L) and Ui 4 ,{T,L) data for the Gaus¬ 
sian ISG in dimension four are shown in Figs. [TTJ [TUJ [TUI 
m The data correspond to = 8192 samples for each 
L. The critical temperature is Tc = 1.80(3) 0 , 0 and 
the finite size critical value for the normalized correla¬ 
tion length ratio [^/L\c = 0.440(5) 0,Q. Scaling against 
the normalized correlation length Fig. [TTl U22{T,L) = 



FIG. 8. (Color on line) Bimodal 3D ISG. Non-self averaging 
parameter U 22 {T, L) against temperature T. L = 4, 6, 8, 10, 
12, 16, 20, 24 (blue, red, black, pink, purple, cyan, green, 
brown) right to left. L = 16, 20, 24 from 0]. Vertical line Tc- 
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FIG. 9. (Color on line) Bimodal 3D ISG. Non-self-averaging 
parameter Uz-i{T,L) against temperature T. L = 4, 6, 8, 10, 
12 (blue, red, black, pink, purple) right to left. Vertical line 
Tc. 

{K 4 ^{T, L)/L)'^ again following the Aharony-Harris law, 
with K 4 = 2.7(1) and peaks located at ^{T,L)/L = 
0.43(2) so very close to ^ITc,L)/L. 

Data obtained for the 4D bimodal ISG (not shown) 
follow a very similar pattern with the same peak height. 
The locations of the C/„„(T, L) peaks are almost inde¬ 
pendent of L. This was noted for U 22 .iT, L) in the 4D 
Gaussian ISG in Ref. [ll|, and in Ref. [l^ for a bond- 
diluted bimodal model; it follows from the proximity of 
the peak ^(T, L) jL and critical ^(Tc, L)/L values. For the 
bond-diluted bimodal model, the peak height is Ri 0.205 
[l 3 l |. Because of the quasi-L-independence, the peak lo¬ 
cation extrapolated to infinite size provides an estimate 
for Tc which is limited in precision only by the statistical 
uncertainties. 

The Gaussian t/„ri(T, L) peak heights become indepen¬ 
dent of L to within the statistical errors after weak fi- 
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FIG. 10. (Color on line) Bimodal 3D ISG. Non-self-averaging 
parameter Ua 4 {T,L) against temperature T. L = 4, 6, 8, 10, 
12 (blue, red, black, pink, purple) right to left. Vertical line 
Tc 



FIG. 11. (Color on line) Gaussian 4D ISG. Non-self-averaging 
parameter 1/22 (F, L) against the normalized correlation length 
^/LT, L)/L. L = 4, 6, 8, 10, 12 (blue, red, black, pink, green). 
The straight line has slope 4. 


nite size effects for small i. Figs. [TH [1311131 The peak 
height values t/ 22 (T,L)max = 0.210(5), [/ 33 (T,L)max = 
0.400(5), C/ 44 (T, F)niax = 0.63(2) are the same as those 
in dimensions 2 and 3 to within the statistical precision. 
The stability of the Unn{T, L) peak heights as L is varied 
turns out to be a useful empirical criterion for the quality 
of equilibration. 

Alternatively, considering the ?7„„(r, L) as dimension¬ 
less variables, the intersections of the curves for hxed L 
should also give a criterion for estimating Tc, but the 
statistical fluctuations and corrections to scaling affect 
the intersections much more drastically than they do the 
peak location, which means that this is an imprecise cri¬ 
terion in the 4D case as noted by 
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FIG. 12. (Color on line) Gaussian 4D ISG. Non-self-averaging 
parameter (722 (F,F) against the temperature T. L = 4, 6, 8, 
10, 12 (blue, red, black, pink, green) from right to left. The 
horizontal line is an extrapolation to criticality at T = Tc 
(vertical line). 
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FIG. 13. (Color on line) Gaussian 4D ISG. Non-self-averaging 
parameter (733 (r,T) against the temperature T.L — 4, 6 , 8 , 
10, 12 (blue, red, black, pink, green) from right to left. The 
horizontal line is an extrapolation to criticality at T = Tc 
(vertical line). 


VI. DIMENSION 5 

U 22 {T,L),U 33 {T,L) and Ui 4 ,{T,L) data for the Gaus¬ 
sian ISG in dimension five are shown in Figs. [T31 UHl HZl 
[TSl The data correspond to 4096 samples for each L. 
The critical temperature is Tc = 2.390(5) and the finite 
size critical value for the normalized correlation length 
ratio [C/T]c ~ 0.45 M- We are not aware of other 
comparable simulation measurements in this dimension. 
Data obtained for the 5D bimodal ISG (not shown) are 
very similar. The (7„„(T, L) peak heights become inde¬ 
pendent of L to within the statistical errors after weak 
finite size effects for small L. The peak height val¬ 
ues U 22 {T,LU.. = 0.215(5), [/ 33 (T,L)^ax = 0.405(5), 
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FIG. 14. (Color on line) Gaussian 4D ISG. Non-self-averaging 
parameter 1 / 44 (T, L) against the temperature T . L = 4, 6 , 8 , 
10, 12 (blue, red, black, pink, green) from right to left. The 
horizontal line is an extrapolation to criticality at T = Tc 
(vertical line). 



FIG. 15. (Color on line) Gaussian 5D ISG. Non-self-averaging 
parameter U 22 {T, L) against the normalized correlation length 
^{T, L)/L. L = 4, 6 , 8 , 10 (blue, red, black, pink). The 
straight line has slope 5. 


C /44 (T,i)max = 0.64(2) are again practically the same as 
those in dimensions 2, 3 and 4 to within the statistical 
precision. 

When scaled against the correlation length ratio, in the 
paramagnetic range U 22 {T, L) = L)/L]^ follow¬ 

ing the Aharony-Harris law Q, with = 2.5(1). The 
peak locations correspond to ^{T,L)/L ft! 0.50. As this 
ratio is greater than [^/A]c the locations of the Unn(T, L) 
peaks are at temperatures below Tc and the peak tem¬ 
peratures move upwards towards Tc with increasing L. 
The peak location extrapolated to infinite size provides 
an estimate for Tc which is again limited by the statis¬ 
tical precision but which provides a useful independent 
check on the value of the ordering temperature. 
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FIG. 16. (Color on line) Gaussian 5D ISG. Non-self-averaging 
parameter U 22 (T, L) against the temperatureT. L = 4, 5, 6 , 7, 
8 , 9, 10 (blue, green, red, olive, black, orange, pink) from right 
to left on the right. The horizontal line is an extrapolation to 
criticality at T = Tc. (vertical line) 
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FIG. 17. (Color on line) Gaussian 5D ISG. Non-self-averaging 
parameter Ud,i(T, L) against the temperature T.L = 4, 5, 6 , 7, 
8 , 9, 10 (blue, green, red, olive, black, orange, pink) from right 
to left on the right. The horizontal line is an extrapolation to 
criticality at T = Tc. (vertical line) 


VII. DIMENSION 7 

By this dimension, N the number of spins per sample 
has become very large, {N = 823, 543 for L = 7), which 
imposes practical limits on the sizes and numbers of sam¬ 
ples in the simulations. The simulations were carried out 
for L = 3 to 7 with 512 samples at each L. 

The dimension 7 bimodal ISG has an ordering tem¬ 
perature Tc = 3.39(1) estimated using the standard 
Binder cumulant crossing point technique [l^ in agree¬ 
ment with the high tem per ature series expansion (HTSE) 
estimates Tc = 3.37(2) [H and Tc = 3.384(15) [3 . (Cu¬ 
riously the HTSE value given in Ref. corresponds 
to Tc = 3.459. We suspect a typographical error). As 
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FIG. 18. (Color on line) Gaussian 5D ISG. Non-self-averaging 
parameter Uaa{T,L) against the temperature T. L = 4, 5, 6, 
7, 8, 9, 10 (blue, green, red, olive, black, orange, pink) from 
right to left on the right. The horizontal line is an extrapola¬ 
tion to criticality at T = Tc. (vertical line) 

this dimension is above the upper critical dimension 
d = 6 , the critical exponents 7 = 1 and zz = 1/2 are 
known exactly. In this case in the paramagnetic regime 
C/ 22 (7’, L) = L) /Tv)®, with an exponent which ap¬ 

pears to be « 6 rather than 7, Fig. [121 This could arise 
from the breakdown of the relations between scaling ex¬ 
ponents above the ucd. Because of the limited number 
of samples and the small values of L at this dimension, 
this estimate is not very precise. 

In the plot of U 22 {T,L) against T, Fig. [201 the L- 
independent critical finite size crossing point value is 
U 22 {Tc) ~ 0.15, and the [C/ 22 (T, T)]max peak heights are 
independent of L and equal to « 0.21 to within the statis¬ 
tics, as for the other dimensions. The maxima locations 
move towards Tc from within the ordered regime. This 
behavior is very similar to that observed in the mean held 
ISG SK model [II [H. 

The higher order U^^{T,L) and C/ 44 (T, T), Fig.[2T]and 
Fig. 1221 follow much the same pattern, with peak maxima 
of 0.39(1) and 0.61(2) respectively, again equal to the 
values for the other dimensions to within the statistics. 


VIII. THE GAUGE GLASS 

The Gauge glass (GG) is a vector spin glass which does 
not support chiral ordering. The GG in dimension 3 has 
a critical temperature Tc = 0.47(1) [20l - [^ . The non¬ 
self-averaging parameter C/ 22 (T, T) scales with ^(T, T)/T 
and shows a maximum peak height independent of L 
and a peak position Tmax(T) near ^(T, T)/L = 0.35. The 
paramagnetic regime data [ 2 ]| appear by inspection to 
be compatible with the Aharony-Harris rule C/ 22 (T, T) ~ 
[^(T,T)/T]® although the published data are not pre¬ 
sented in this way. As the critical correlation length ratio 
is [^(T,T)/T]c = 0.54 [ 2 l|, the C/22(T,T) peak tempera- 



FIG. 19. (Color on line) Bimodal 7D ISG. Non-self-averaging 
parameter C /22 (T, L) against the normalized correlation length 
^(T,L)/L. L = 4, 5, 6, 7 (blue, green, red, olive). The 
straight line has slope 6. 
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FIG. 20. (Golor on line) Bimodal 7D ISG. Non-self averaging 
parameter C /22 (T,L) against temperature T. L — 3, 4, 5, 6, 7 
(cyan, blue, green, red, olive) right to left. Vertical line Tc- 
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FIG. 21. (Golor on line) Bimodal 7D ISG. Non-self averaging 
parameter C/ 33 (T, T) against temperature T. L — 3 , 4, 5, 6, 7 
(cyan, blue, green, red, olive) right to left. Vertical line Tc. 
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FIG. 22. (Color on line) Bimodal 7D ISG. Non-self averaging 
parameter Uaa{T, L) against temperature T. L = 3, 4, 5, 6, 7 
(cyan, blue, green, red, olive) right to left. Vertical line T^. 

ture location moves downwards with L and tends towards 
Tc- The 3D vector spin glass GG U 22 {T,L) thus follows 
basically the same rules as followed by 1/22 (T,L) in the 
Ising spin glass in 3D, except that the GG peak maximum 
is ~ 0.10 instead of 0.205. Data on GGs in dimensions 
2, 3 and 4 from measurements which were not designed 
to estimate the non-self-averaging parameter are con¬ 
sistent with U 22 {T,L) peak values near 0.10 in each di¬ 
mension. We can speculate that this family of spin glass 
models also has its characteristic dimension-independent 
value of the non-self-averaging peak height. 

IX. HEISENBERG AND XY SPIN GLASSES 

Numerical measurements on Heisenberg spin glasses 
(HSGs) are of particular importance because the canon¬ 
ical experimental spin glass dilute alloys (AuFe, CuMn, 
AgMn) are all Heisenberg systems, so it should be pos¬ 
sible to understand the ordering mechanism in ” real life” 
spin glasses on the basis of numerical data on Heisen¬ 
berg models. We have no new data to report on these 
models but it is of interest to consider published non-self¬ 
averaging data in the light of the ISG results. 

Both Heisenberg and XY spin glasses can support chi¬ 
ral glass order as well as spin glass order, and for many 
years there have been two conflicting interpretations of 
the numerical data on the ordering transitions in these 
models in dimension 3. According to the first interpre¬ 
tation, the ordering is spin-spin interaction driven; basi¬ 
cally the ordering process is much the same as in ISGs, 
and the chiral order follows on as a geometrically neces¬ 
sary consequence of the onset of spin order, without the 
chiral interactio ns playi ng any significant role in the spin 
glass transition [23l - l27l |. The alternative interpretation 
is that the driving role in 3D HSG or XYSG ordering is 
played by the chirality, so that there is first a chiral or¬ 
der onset followed at a lower temperature by spin order¬ 


ing transition [28l - l3lj |. (Similar disagreements concerning 
fully frustrated 2D XY models were resolved definitively 
in favor of a distinct chiral order transition just above a 
spin order transition [3^ . Hi). The arguments of both 
schools to support their respective interpretations in the 
3D HSG and XYSG models have been essentially based 
on analyses of the data for the crossing points of the 
dimensionless normalized spin and chiral (s and c) corre¬ 
lation lengths ^s{T, L)/L and ^c{T^ L)/L. The numerical 
simulations in the spin glasses are even more demanding 
than in the fully frustrated models, and because of intrin¬ 
sic finite size corrections and the need to reach strict equi¬ 
libration at each L, extrapolations to infinite L in order 
to estimate the ThL crossing point locations are delicate. 
As simulations were extended to larger sizes in successive 
Gaussian HSG and XYSG measurements interpreted on 
the spin-driven ordering scenario, the joint spin/chiral 
crossover temperature was estimated to be Tc{HSG) ~ 
0.160 dlL Tc(ff5'G) Ri 0.145 with a KTB-like criti¬ 
cal line |24j . marginal but very similar spin and chiral 
behavior (XYSG and HSGjflH, H^, and most recently 
Tc{HSG) w 0.120 [ 13 . No non-self-averaging results 
were reported. From detailed 3D bimodal and Gaussian 
HSG and 3D Gaussian XYSG measurements, the two 
separate transition temperatures on the chiral-driven or¬ 
dering scenario are estimated to be (bimodal HSG) [^ . 
Tc(c) = 0.194(5) and Tc(s) < 0.15, (Gaussian HSG) 
^c) = 0.143(3) and Tc(s) = 0.125(-f0.006/ - 0.012) 
m, and (XYSG) Tc(c) = 0.308(5) and T^{s) = 0.274(3) 
3l| . Non-self-averaging data were shown in each case. 

In the light of the ISG results reported above, it would 
appear that in Heisenberg and XY spin glasses non-self¬ 
averaging could provide an independent primary numer¬ 
ical criterion for spin and/or chiral ordering much less 
sensitive to hnite size effects and to strict equilibration 
(as already suggested in Ref. [l^). On the first (spin- 
driven ordering) scenario one would expect the spin non¬ 
self-averaging parameter U 22 s{T,L) to follow much the 
same rules as for the ISG or the GG chiral-free vec¬ 
tor spin glass cases discussed above, with a peak loca¬ 
tion moving towards an upper spin ordering temperature 
Tc(s) as L increases, and a regular behavior reflecting 
U 22 s{T,L) ~ {^s{T,L)/L)^ in the paramagnetic regime 
above Tc(s). On this interpretation the chiral C/ 22 c(T, L) 
would be weaker than the U 22 s{T, L); if a U 22 c{T, L) peak 
exists, it would be located at a temperature below or pos¬ 
sibly at the ThL U 22 s{T,L) peak. 

On the second (chiral-driven order) scenario, it would 
be the chiral U 22 c{T, L) which would show a peak first, 
with a peak location tending towards the (upper) chi¬ 
ral ordering temperature Tc(c) as L increases. In the 
paramagnetic regime one would expect a regular behav¬ 
ior of the chiral non-self-ordering 1722c (T, L) with increas¬ 
ing L, governed by U 22 c{T,L) ^ {^c{T, L)/L)^. On this 
scenario the spin U 22 s{T,L) would then show a peak 
location somewhere below the chiral 1722c (T, L) peak, 
with a location tending towards an ordering temperature 
Tc(s) below Tc{c), together with a paramagnetic regime 
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U 22 s{T, L) behavior behaving irregularly at least at small 
L because the paramagnetic spin ordering is perturbed 
by the dominant onset of chiral order. 

Very informative non-self-averaging data have been 
published on the 3D HSG with bimodal interactions [^ . 
on the 3D HSG with Gaussian interactions [s^, and on 
the 3D Gaussian XYSG [3l|. In each case the pattern is 
the same : 

- there is a strong U 22 c{T,L) peak at an almost L- 

independent temperature T « 0.19, T Ri 0.145, T « 0.31 
respectively, so in each case close to the Trie) value es¬ 
timated independently from other criteria In 

the paramagnetic regime there is a regular narrowing 
in temperature of the U 22 c{T, L) peak with increasing L 
which appears compatible with the Aharony-Harris law 
U 22 c{T,L) ^ {^c{T, L)/L)^ though the data are not pre¬ 
sented in this form. 

- in each case, the spin U 22 s{T, L) peak is either not vis¬ 
ible (HSGs) or is marginal (XYSG) down to the lowest 
temperature at which non-self-averaging measurements 
were made, T Ri 0.145,T r; 0.11 to 0.13 depending on 
L, and T r 0.24 to 0.275 depending on L in the three 
cases. Over the whole temperature range U 22 s{T, L) is 
always weaker than U 22 c{T, L), and U 22 siT,L) has ir¬ 
regular behavior as a function of L in the paramagnetic 
temperature regime at and above the U 22 c{T, L) peak. 

Thus the non-self-averaging data U 22 c{T,L) and 
U 22 s{T,L) in the three models [2^-[3lj can be seen by 
inspection to be clearly incompatible with the behav¬ 
ior expected on the spin-driven ordering scenarios 

, and fully compatible with the chiral-driven ordering 
interpretation [ 2 § - [^ . 

X. CONCLUSION 

The non-self-averaging data on ISGs in all dimen¬ 
sions show a remarkable regularity. In each dimen¬ 
sion there is a peak as a function of temperature 


in the standard non-self-averaging parameter 1/22 (T,L) 
and in the higher order parameters 1/33 (r,L) and 
1/44 (r,L) whose values are L-independent after weak 
small size effects; the peak values C/ 22 (T,L)max ~ 0 . 21 , 
D 33 (r, L)max ~ 0.40, C/ 44 (T,T)max ~ 0.62 are indepen¬ 
dent of dimension to within the statistics for dimensions 
2,3,4, 5 and 7. In the paramagnetic regime above the 
peak the Aharony-Harris renormalization group law [l| 
Unn{T,L) = {Kd^{T,L)/LY is obeyed, with Kd{U 22 ) ~ 
2.6 for all dimensions. Both of these empirical observa¬ 
tions can be classed as “hyperuniversal behavior”. 

Published [2l[ and unpublished data on the Gauge 
Glass, a vector spin glass which does not support chi¬ 
rality, suggest that non-self-averaging rules analogous to 
those that hold in the ISGs appear to apply but with a 
different characteristic peak height t/ 22 (T, L)max ~ 0.10. 

XY and Heisenberg spin glasses can support chiral or¬ 
dering as well as spin ordering. In the light of the non¬ 
self-averaging behavior reported above for the ISG mod¬ 
els, it is clear that the published spin and chiral non¬ 
self-averaging data [29l - l3l]| in 3D Heisenberg and XY 
models are incompatible with a spin-driven ordering sce¬ 
nario [23 I - [t^ but strongly support the alternative con¬ 
clusion that the spin glass ordering in these models is 
chiral-driven rather than spin-driven, on the Kawamura 
scenario [2^. An important implication is that order in 
the canonical experimental Heisenberg spin glasses is also 
chirality driven. 
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